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Heavy fermion superconductivity: X-boson treatment 
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We use the X-boson method to study the heavy fermion superconductivity phase employing an 
extension of the periodic Anderson model in the U — oo limit with a nearest neighbor attractive 
interaction between the localized /-electrons. We show that higher values of the hybridization 
parameter implies in lower values associated to the maximum of the superconducting critical tem- 
perature T c , indicating that the real charge transfer between bands tends to destabilizes the Cooper 
pairing. Moreover, we show that the superconductivity is constrained to the vicinity of a range of 
densities where the /-band density of states at the Fermi level Pf(fi) have sufficiently high values 
and it was found both for configurations where the system presents intermediate valence and heavy 
fermion behavior, as experimentally observed. Finally, as the total occupation is raised and for 
larges hybridization parameter V there is an insulator-superconductor transition, which is related 
to the existence of a hybridization gap that cannot be accessed by the slave boson method, because 
when the chemical potential lies in or above the gap the system suffers the second order phase 
transition characteristic of that method. 

PACS numbers: 74.70.Tx, 74.20.Fg, 74.25. Dw, 75.30. Mb 



I. INTRODUCTION 

In the present paper we study the heavy fermion su- 
perconductivity employing an extension of the periodic 
Anderson model (PAM) in the U = oo limit taking into 
account the effect of a nearest neighbor attractive in- 
teraction between / electrons The model was 
primarily designed to study the heavy fermions systems, 
but it could fit in a description of the high temperature 
superconductors compounds (HTSC) as well, if regarded 
as a simplified extension of the Emery model 0, |S[ . 

Heavy fermion materials present a great variety of 
ground states: antiferromagnetic (AF) (UAgCu4, UCU7), 
superconductors (CeCu2Si2, UPt 3 ), Fermi liquids (FL) 
(CeCu6, CeAW) and Kondo insulators (KI) (Ce3Bi4Pi 3 , 
YbBi 2 ) H, H EJ • An uniform high temperature Curie- 
like magnetic susceptibility, a common feature to these 
compounds, is related to the fact that they are composed 
of elements with the incomplete /-shells like Ce and U. 
As the temperature decreases up to a certain range, the 
system presents a temperature independent uniform sus- 
ceptibility (Pauli susceptibility) signaling the annihila- 
tion or binding of the magnetic moments of the localized 
/-states, resembling the single-impurity Kondo problem 
|10|. A consistent description of the overall properties 
of the heavy fermions is achieved by the competition be- 
tween the Kondo effect, dealing with the annihilation or 
binding of the localized magnetic moments, and the Rud- 
erman Kittel Kasuya Yosida (RKKY) interaction, which 
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favors the appearance of a magnetic ground state; thus, 
the Hamiltonian that describes the basic physics of the 
system may be a regular lattice of /-moments which in- 
teract with an electron gas and with themselves through 
these electrons as proposed by the PAM. As concerns the 
heavy fermion superconductivity, it is an experimental 
fact that the large specific heat jump at T c , of the order 
of the normal phase specific heat, indicates that pair- 
ing takes place among the /-electrons and the coherent- 
length must be much shorter than the typical values 
for conventional superconductors @. Furthermore, the 
narrow bandwidth of the /-electrons suggests a strong 
on-site Coulomb repulsion which precludes on-site pair- 
ing and may lead to non-s-wave or unconventional pair- 
ing. New results have also revealed the possibility that 
magnetic (mainl y A F) phase may coexists with the su- 
perconductivity [llj and that some heavy fermion com- 
pounds display pseudogap behavior at the normal phase 
[12fl which is a common feature of the high critical tem- 
perature superconductors |l3j. 

These rich phenomena motivate us to study the su- 
perconducting phase in context of the PAM with the re- 
cently developed X-boson technique, that was applied to 
the PAM in the limit U -> 00 [HEaE^- This a PP roach 
was partially inspired in the mean-field approximation of 
Coleman's "slave boson" method 0, 0] and produces 
satisfactory results, because they are very close to those 
obtained by the slave boson in the Kondo limit at low 
temperatures, while it recovers those of the chain cumu- 
lant approximation (CHA) [H I |2C | . which is deriv- 
able [21], [2^ and is a very good approximation at high 
temperature T. The unphysical second order phase tran- 
sitions that appear in the slave boson approach when 
fj, » Ef.ja = Ef at low T, fi being the total chemi- 
cal potential and Et the energy level of the /-electrons, 
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and also for all parameters at intermediate temperatures 
T, are then eliminated by the X-boson treatment. Cole- 
man 0| has observed that these phase transitions are 
artifacts of the theory, and the advantage of the present 
treatment is that those spurious phase transitions do not 
occur. Although in the X-boson there is no relevant spa- 
tial dependence in the self-energies of the CHA, they do 
retain some local time dependence. In the mean-field 
slave boson the corresponding self-energy vanishes, show- 
ing that all the local time dependence is completely lost 
in the self-energies of that method. The spurious phase 
transition observed by the slave-boson method are absent 
from the X-boson description, and this result indicates 
that the time dynamics retained in the CHA is able to 
suppress those transitions [lfij . 

In the present paper we adopt the schematic classifica- 
tion proposed by Varma |23j and recently reintroduced by 
Steglich et al [(J 13 *° describe the exemplary Ce-based 
compounds. It is given in terms of the dimcnsionlcss cou- 
pling constant for the exchange between the local / spin 
and the conduction-electron spins, g = Np\J\, where Np 
is the conduction-band density of states at the Fermi en- 
ergy and J is connected to the parameters of the PAM 
via a Schrieffer- Wolff transformation |25|. J = 2V 2 /Ef 
when U — > oo. Therefore, within the X-boson technique, 

g = p D (jj)2V 2 /\E f \, (1) 

where p c {p) is the conduction density of states at the 
chemical potential p, V is the hybridization parameter 
and Ef is the renormalized Ef introduced by the X-boson 
method, as shall be explained below, and the behavior of 
these compounds can be qualitatively driven through this 
parameter: when g > 1, the compound under considera- 
tion presents an intermediate valence (IV) behavior while 
for g < 1 it presents heavy fermion Kondo regime (HF). 
There exists a critical value g c at which the Kondo and 
the RKKY interactions have the same strength and non 
Fermi-liquid (NFL) effects have been postulated for sys- 
tems with g = g c . For g c < g < 1, the magnetic local 
moments do not exist at very low temperatures and the 
system presents a Fermi liquid behavior while for g < g c 
we have the local magnetic moment regime. We observe 
that we use the parameter g to classify the regimes of the 
PAM in a very qualitative way. Finally, as concerns to the 
X-boson chain approach, its present form only includes 
hybridization effects to second order in V, therefore the 
RKKY effects are not taken into account and we cannot 
discuss non Fermi liquid behavior at the present stage of 
the work. 

The paper is organized as follow: In Section[I]we intro- 
duce the main ideas of heavy fermion superconductivity, 
in Section [H] we make a brief revision of the X-boson 
method developed earlier ^5|, in Section ITTT1 we present 
the model, in Section IIVI we calculate in a mean field 
approximation the Gorkov's anomalous function written 
in Zubarev's notation and the gap equation, in Section 
[V]we obtain numerically the superconducting phase dia- 



gram for several hybridization parameters and we discuss 
the model results; finally in the Section lVIl we present the 
conclusions and the future directions of the work. 



II. THE X-BOSON METHOD 

The periodic Anderson model (PAM) in the limit of 
infinite Coulomb repulsion U — > oo is given by 

H = ek, CT c ki0 .Ck,<7 + ^2 E f,i° X i,™ 

k,fT jjCr 

+ £ (r / ,, k .v;,„ T r k .. + \ :, k r k , r x !J:1 ) , (2) 

where the first term is the kinetic energy of the conduc- 
tion electrons (c-electrons) , the second term describes in- 
dependent localized electrons (/-electrons) of energy Ef, 
and the last term is the hybridization Hamiltonian giv- 
ing the interaction between the c-electrons and the /- 
electrons. 

We employ the Hubbard operators |2y] to project the 
double occupation state | j, 2) with two local electrons 
out from the space of local states at site j, as the X Hub- 
bard operators do not satisfy the usual commutation re- 
lations, the diagrammatic methods based on Wick's the- 
orem are not applicable, one has to use product rules 
instead: 

Xj.ab-Xj.cd — Sb^cXj^ad- (3) 

The identity decomposition in the reduced space of local 
states at a site j is then 

Xj.00 + A^crer + Aj : ctct = Iji (4) 

where a = — <r, and the three Xj^ aa are the projectors 
into | j, a). Because of the translational invariance, the 
occupation numbers rij, a = < Xj <aa > satisfy Tij >a =n a 
(independent of j), and from Eq. (QJ we obtain the "com- 
pleteness" relation 

no + n a + n w = 1 . (5) 
In Coleman's "slave boson" method 

[HEi, the Hub- 
bard X operators are written as a product of ordinary 
bosons and fermions: Aj 00 — > &T bj, Xj oa — > bjfj jCn 
Xj,ao — > fj~o-bj, and a condition, that is equivalent 
to Eq. JSJ, is imposed to avoid states with two elec- 
trons at each site j. In the mean field approximation 
b~l b~l >= yfz and the method of Lagrangian multi- 
pliers is employed to minimize the free energy subject to 
that condition. The problem is then reduced to an uncor- 
related Anderson lattice with renormalized hybridization 
V — * y/zV and / level ej — » e/ + A, and the conservation 
of probability in the space of local states is automatically 
satisfied because they are described by Fermi operators. 

The approximate GF obtained by the cumulant expan- 
sion |2jj|30j do not usually conserve probability (i.e. they 
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do not satisfy Eq. (0), and the procedure we adopt to 
recover this property in the X-boson method is to intro- 
duce 

R =< Xj, 00 >=< bjbj >, (6) 

as variational parameter, and to modify the approximate 
GF so that it minimizes an adequate thermodynamic po- 
tential while being forced to satisfy Eq. (0. To this pur- 
pose we add to Eq. (|2J) the product of each Eq. Q into 
a Lagrange multiplier Aj , and employ this new Hamilto- 
nian to generate the functional that shall be minimized 
by employing Lagrange's method. To simplify the calcu- 
lations we use a constant hybridization V, as well as site 
independent local energies Efj^ a = Ef^ a and Lagrange 
parameters Aj — A. We then obtain a new Hamiltonian 
with the same form of Eq. J3J : 

+N S A{R-1), (7) 
but with renormalized localized energies 

Ef, a = e fta + A. (8) 

The parameter 

R = 1 - ]T < X aa > (9) 

cr 

is now varied independently to minimize the thermody- 
namic potential, choosing A so that Eq. (J3J be satis- 
fied. While at this stage the electrons in the slave bo- 
son Hamiltonian have lost all the correlations, the Eq. 
(0) is still in the projected space and it is not neces- 
sary to force the correlations with an extra condition. 
On the other hand we do not have an exact solution for 
this new problem, and we then consider the most simple 
approximation obtained within the cumulant formalism, 
the Chain approximation (CHA) 0,[2(j. The need of 
minimizing a thermodynamic potential arises because the 
completeness relation is not automatically satisfied for 
approximate cumulant solutions, and although the two 
procedures are formally very similar, they have a rather 
different meaning. 

III. THE SUPERCONDUCTIVITY MODEL 

In order to have a superconducting state, we add to 
the usual PAM given by Eq. © , an effective interaction 
Hw among the heavy /-electrons, 

^ w = 2 WijX iiaa Xj^i a i , (10) 

<i,J>,<T,cr' 



where 

Xi t aa = Xj 0(y Xifi a . (11) 

The term Hyy describes an effective attraction between 
two neighboring / sites (Wij < 0), which is responsible 
for the heavy fermion superconductivity in the model. 
We only consider the superconductivity arising from /- 
electron pairing. As the conduction electrons exist at 
the Fermi surface the possibility of formation of c-c or 
f-c Cooper pairs cannot be excluded but we do not 
consider these pairs here. Indeed, there are experimen- 
tal evidence that pairing occurs between the heavy /- 
electrons. In this framework the correlated /-electrons 
hybridize with the conduction band, and interact via a 
non-retarded nearest neighbor attraction. The interac- 
tion is proposed on a phenomenological basis and its mi- 
croscopic origin is not investigated. The same approach 
was also used by Romano et. al. Q, Zielihski and Mat- 
lak 0, and Araujo et. al. 0. Moreover, Tachiki and 
Maekawa [2i| studied the superconductivity in the Pe- 
riodic Anderson Model with a small dispersion of the 
/-band in the heavy fermion state and they have con- 
cluded that pairing between /-electrons are responsible 
for superconductivity, rather than between the conduc- 
tion electrons. 

It is more interesting to work with the Fourier trans- 
form of the Hubbard operators 

X i:0 a = —7= ^ e lk Rl X ki0c r, (12) 



xi„ = -^j2 e - ik - Rl ' x U, (13) 

where N is the number of the lattice sites. Considering 
Cooper pairing only in the singlet channel, the model 
defined by Eqs. J7J and (|10fl in the momentum space 
can be written as 

H = e k,<r c k, a c k,cr + ^ £'/A"k,<j<j 

k,(T k, CT 

+ E y ( X k,o CT Ck, ff + c k , ff AW) 

k,cr 
k,k' 

+N S A(R-1), (14) 

where 

bi = Xl, 0a Xi ktOff (15) 

and 

W klk > =w£V( k " k H (16) 

s 

where the summation over 5 runs over the nearest neigh- 
bors and we consider the hybridization constant V = 
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The superconducting term in the right hand side of the 
Eq. H14|) is responsible for the Cooper pair formation. 
The parameter W is negative, where Wkk' = Wrfcrfc' 
and f?k is assigned according to the symmetry of the su- 
perconducting order parameter Q . 

For most of the superconducting materials the charge 
carriers can couple in the s, p, d, etc. channels. For heavy 
fermions, although, both the crystalline anisotropy and 
spin-orbit interaction are important |3l| . and the order 
parameter should be written in terms of a complete set of 
basis-function multiplets for the appropriate symmetries. 
However, as a first study of the effect of the hybridization 
over the superconducting phase diagram and how the 
hybridization affects T c , we only consider the case of an 
isotropic s-wave superconducting gap in our numerical 
calculations. 



IV. THE GREEN'S FUNCTIONS 

Since our model Hamiltonian, Eq. 114(1. does not de- 
pend explicitly on the time, the corresponding Green's 
functions GF are only functions of the difference t — t' = 
t. In the aforementioned Zubarev's notation |3^ | the 
GF can be written as G^ c (k,r) = 9{t)(c k<7 (t)cI !<t (0)) - 
6(— t)(c] c CT (0)ck, a (t)) for this fermionic system, where 
the operator c k (t) is in the Heisenberg representation. 
Formally, besides the (r, t) representation of the GF, the 
model admits the (r, u>) representation defined by the 
Fourier transform G cc CT (r,u;). Hence the (k, a;) represen- 
tation is defined by the discrete Fourier transform 



G cc .(r,c) = l^^ c (k, W )e 4k - r , 



N 

k 

which could be simply denoted by 



(17) 



(18) 



Analogously, the Gorkov's anomalous function is defined 
by 

FjfJKu;) =« Xt , 0CT ;Xt k0 _ > w . (19) 
In a similar way one defines 

Q S J (-k, u) =« X- K0W ; Xl k _ » w , (20) 



and 



g c J{-^ u) =« c_ k ,a, Xl k » w , (21) 



FlfA*> w ) s « X -Ko* »- • ( 22 ) 



^From the above definitions, one can show that 
a (k, oj) satisfies the equation of motion 



where Zk.cr = X^ 0a ,H and 77 is given by Eq. ijTI)) . 

We employ the equation of motion method and we obtain 
a chain of equations that shall be broken by a mean- 
field approximation and further calculations furnish the 
quasi-particles spectrum. Therefore, the final system of 
equations to be solved is 



(iuj n - ek)0^(-k,w n ) 
(iuj n + ek)Tl }(7 {\L,uj n ) 



(ILU 



V^'C-k.Wn), 

-yjFt /CT (k, w „), 
-yA^^k,^) 



k, uv, 



(24) 
(25) 



(26) 



(iw„ (-k,w B ) = VD^^(-k,w n ) 

+A fe £V^) /iff (k,w n ) 
+Af, (27) 
where A& is the superconducting gap defined by 

Ak = ]rWW(&l> (28) 



kk' 



and the quantity D a is given by 
D a = (X 00 + X 

a 



(29) 



In the X-boson method employed here all the correla- 
tions are included in the quantity D a and, indeed, mak- 
ing D a — 1 the above system of equations with the map- 
ping V <-> \fzV and W <-*■ z 2 W is identical to the slave 
boson's result given in Appendix^ similar to those ob- 
tained by Araujo et. al. 0- Further, in the absence 
of the hybridization term the /-band and the conduc- 
tion band decouples and solving the remaining system of 
equations for the /-band one formally recover the BCS 
result though with a distinct physical meaning, here the 
energy Ef assumes a constant value and the strong cor- 
relations presented are taken into consideration via the 
parameter D a . 

Considering the paramagnetic case (D a = D w ) and 
solving Eqs. 1|24 )I - 1|27JI for the anomalous GF one find 



i=i 



where 



Bo 



B 4 = A] 



(^ 2 -e k 2 ) 
' 2bJ2 {uJ2 2 — ^4 2 ) 

Da 2 (^4 2 - 6k 2 ) 
' 2W4 (<x>4 2 — Ll>2 2 ) ' 



(30) 



(31) 



(32) 



with B\ = —B2 and B3 = —B4. The four poles repre- 
senting the above quasi-particle (QP) excitations given 
by 



(23) 



W2 



U>4 



4/3, 



4/3, 



(33) 
(34) 



5 



u>x — —u>2 and W3 = W4, where 



a = A k + 2V 2 D a + e k 2 + Ep , 



(35) 



and 



= V"D a l - 2V z D a e k E F + A k 'e k z + e k z E F . (36) 

Now we shall treat the hybridization and the supercon- 
ducting terms (Hy and Hw respectively) in the model 
Hamiltonian given by Eq. I|14|) as external perturbations 
and calculate the thermodynamic potential £1 through 
the A parameter integration 33]. One should minimize 
the free energy in order to find the Lagrange multiplier 
in the former Hamiltonian. But first we will focus on the 
contribution of Hw only, which is the only term respon- 
sible for the Cooper pair formation. Our discussion is 
analogous to imposing V = in our model. Under this 
constraint the difference of thermodynamic potential in 
the superconductor and normal states tt s — tt n is 



n s -n n = f d\(H{ 

Jo 



where 



Hi 



w 



k,k' 



IV / 



'bib* 



k'O k 



(37) 



(38) 



By a variable transformation the above difference defined 



by Eq. ^3 is rewritten as $™ {dW /W'){H' W ) . In a 
Hartree-Fock approximation the average (H^) can be 



"factorized" and 5^ k , r/k^k'^k' 
Hence Eq. l|37|l becomes 



A k /VT according to Eq. 



k 



W 



dW 



A k <fe[; 



(39) 



For an isotropic s-wave superconducting gap in the ab- 
sence of an external magnetic field A k = A = A* and 
r/ k = 1. Also notice that W~ 2 = -dW' 1 /dW and after 
another variable transformation we get that 



dA' A" 



dA' 



(40) 



which is the same functional result given by Fetter and 
Walecka [34[ • Notice that at temperatures higher than or 
equal to the superconducting critical temperature T > T c 
the superconducting order parameter is assumed, as in 
the BCS theory, to be zero. Therefore il s — Q„ = at 
T = T c and henceforth we shall disregard the contribu- 
tion of Hw in order to calculate the Lagrange multiplier 
for the superconducting phase diagram and only consider 
the hybridization term contribution of the model Hamil- 
tonian as an external perturbation. Further, for T > T c 
the system is in the normal state and the GFs yields to 
the previous result obtained in the CHA of the PAM, 
obtained by taking the infinite sum of all diagrams that 



contains only second order cumulants terms. The GFs of 
the CHA have functional form which are close to the un- 
correlated ones (U — 0), but they cannot be reduced to 
them by any change of scale, except for D a = 1, when we 
recover the slave-boson GFs if we put V — + \pzV . Nev- 
ertheless, notice that although the condition that forces 
completeness in the CHA is identical to that employed in 
the slave-boson method to force nt < 1, it has a rather 
different origin, being only a consequence of using a re- 
duced set of diagrams in the perturbative expansion, and 
the departures from completeness are usually very mod- 
erate. In the formalism described here, it is this essential 
difference between the two methods that eliminates the 
spurious phase transition appearing in the slave-boson 
method. Imposing A = at the system of equations 
given by Eqs. pi )) -l'2*7 |) one find that 



Qf f (-k,u n ) \t=t c = 



B 



lUJ n — U) + lOJ n — LO- 

where the poles w+, lo~ are 



(41) 



AV 2 D a 



z 

and the above coefficients are 

e k - w 4 



0J- - £k 



A = D a 
B = D C 

1*)+ — UJ— 

Also notice that 

(iw tt -ek)V- 1 eS / (-k,w») =Q f J{-^n) 
and one get that 

VD„ 



(42) 

(43) 
(44) 



G c J{-k,w n ) 



1 



1 



UJ + lLU n 



Moreover, 
G c *(k,uj n 



Ef 



(iui n - Ef)(iu n - ek) - V 2 D a 



(45) 



(46) 



Since the thermodynamic potential was already calcu- 
lated in detail by one of us and collaborators elsewhere 
|15| treating the hybridization Hamiltonian Hv as an ex- 
ternal perturbation, we shall just briefly outline its so- 
lution below. According to the A parameter integration 



n = Q + N s A{R-l)+ dX(H v {X)) x , (47) 

Jq 

where, regarding Eq. 114*51) . the average (Hv(X))\ is 

(Hy(X)h = iX^D ™^ , (48) 
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and riF is the Fermi function. At the above Eq. (|47|l the 
quantity O is the thermodynamic potential associated 
to the unperturbed part of the Hamiltonian at T = T c 
and iV s A (R — 1) was already introduced to impose the 
completeness relation Q. Minimizing SI with respect to 
R, we get 



nF(w+) — np(toJ) 



(49) 



e k - E f ) + AV 2 D„ 



Finally, the average can be found from Eq. 130(1 
and the self-consistent gap equation is obtained accord- 
ing to the definition in Eq. J5SJ. Hence, at the supercon- 
ducting critical temperature the superconducting gap for 
an isotropic s-wave superconductor satisfies the relation 



i = -wp- 1 Yl ■ 



Di 



n.k 

X 



[<4+u£ (2V 2 
(v 2 -e k E f y 



^ 2 



4 + Ef 



Notice that Eq. (|50jl is identical to 

n , k 



T=T C 



(50) 



(51) 



^From the numerical point of view it is convenient to 
perform first the k summation in the Eq. I|51|) and we 
get that 



1 = -- 



WD„ 



1 



1 



with 



ai 



V 2 DgE f 



(52) 



(53) 



and 



For simplicity we only consider a constant conduction 
density of states in numerical calculations, 

p(e k ) = lp> ~ D <^-»< D , (55) 
rv ' I , otherwise v ' 

and integrating Eq. i|52|) over the above square band 
l)55p. we get that 



1 WD °<t 

1 — o 

Pc 



(56) 



where 



Sm — 



1 
2~D 



arctan 



D-ai 



u 2 n + (al-bj) 
h 

D + a\ 



arctan 



bi 



2L> 



In 



6f + (D + ai) 



(57) 



The self-consistent solution of the Eqs. itSSl 
and l|57|l provide the superconducting phase diagram. 

In the next section we calculate the superconducting 
phase diagram for T c as a function of the total electron 
number N t = Nf + N c , where these occupation numbers 
are calculated for the Green's functions given by Eqs. 
(|41|l . I|45(l and (|46() considering a constant conduction 
density of states as defined in Eq. Q55fl . 

The results are 



LO — E f 



DgV 

oj — Ef 



# C M = 

where 



x log 



log 

2D 6 



VD. 



uj-E 



[cj — U}_ 



(58) 
(59) 
,(60) 



: \UE f + D. ) +W, T V- 
(Ef+D. 



1 Mx. 
2 



± 



with D± = D ± [i. 
Also, 



AD a V 2 



Nr. 



dun F (u)Q[g^(u)] , 



(61) 

(62) 
(63) 

(64) 
(65) 



where n F is the Fermi function. 



V. RESULTS AND DISCUSSION 

In this section we solve the self-consistent system of 
equations Eqs. ©, gHJl, (JSSJl and for T c . All the 
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energies are expressed in units of D, numerical calcula- 
tion are performed making Et = —0.15 and the super- 
conducting interaction parameter is W = —0.10. 

In Fig. n we present the superconducting phase dia- 
gram for T c as a function of N t , where each N t corre- 
sponds to a distinct sample with a certain stoichiomet- 
ric composition and the chemical potential is found self- 
consistently constrained to the furnished N t . The super- 
conducting interaction W is kept constant for every set 
of parameters, since its main effect is simply to raise the 
T c value. 
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FIG. 2: The occupation numbers N c and Nf as a function of 
Nt. Other parameters are the same as in Fig. Q 
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FIG. 1: The superconducting critical temperature T c as a 
function of the total electron number Nt for (a) V = 0.1, (b) 
V = 0.3, (c) V — 0.5 grouped in the three different regimes. 
Ef = —0.15 and the superconducting interaction parameter 
W is W = —0.10, all the energies are expressed in units of D. 

Our results for the superconducting phase diagram 
show that the existence of the superconductivity is con- 
strained to the vicinity of a range of occupations where 
the values of the /-band densities of states p/(/i) at ui = fx 
are sufficiently high. Moreover, high values of usu- 
ally implicate higher values of T c , what suggests that the 
Kondo behavior of the system favors superconductivity 
in the X-boson approach. Nevertheless, the reduction of 
the hybridization parameter V causes a general increas- 
ing on the maximum of the superconducting critical tem- 
perature T c , this is due to the diminution of the charge 
fluctuations between the conduction and the / electron 
states. Indeed, Fig. [21 shows the occupation numbers N c 
and Nf as a function of Nt, for a given total electron 



FIG. 3: The critical temperature T c as a function of the 
hybridization V for Nt ~ 1.4 and Nt = 1.8. Other parameters 
are the same as in Fig. 



number as V increases, part of the electrons that would 
be in the localized /-band, and therefore responsible for 
the Cooper pair formation, should be hybridized with the 
conduction band, showing that the charge fluctuation be- 
tween the conduction and the /-electron states makes the 
Cooper pairs unstable. Similar results were also obtained 
by Romano et al. Q and Sarasua and Continentino |35| . 
Indeed, Fig. shows the behavior of T c as a function 
of the hybridization V in the low and high occupation 
regimes. As V increases, the critical temperatures di- 
minishes and at very small values of the hybridization, 
the critical temperatures tend to a finite value; this be- 
havior is similar to that obtained by Araiijo et. al. Q 
using the slave boson method. Note that for Nt = 1.8 
the superconductivity phase is suppressed for the interval 
V ~ [0.03, 0.06], which is related to the appearance of a 
hybridization gap in the pf(to). 

Regarding the superconducting phase diagram, at low 
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V = 0.1, (d)-(f) V — 0.3, and (g)-(h) V = 0.5 compared to the values of the chemical potential p (vertical short dashed line). 
Other parameters are the same as in Fig. 



total occupation N t , the / level occupancy is also small 
and, since Cooper pairing occurs only between the / 
electrons, T c is null. As the values for Nt are raised, 
T c increases as Pf(p) is increased and presents its maxi- 
mum in the vicinity of the Kondo resonance. In the large 
Nt regime the / level is almost fully occupied with one 
electron per state and the Pf{p) decreases causing the 
suppression of superconductivity, but it should be noted 
that T c vanishes before Nf — > 1. Similar results were 
also obtained by Araujo et al 0, except that for them 
the superconducting critical temperature is constrained 
to T c < Tk, where the Kondo temperature Tk is de- 
fined as the T that makes the slave-boson parameter z 
vanish. Indeed, in the slave-boson method by increasing 
the temperature T or the chemical potential p, the value 
V = \/zV = is presently reached, leading to Nf — > 1 
and the decouple of the two types of electrons, what can 
be interpreted as a change of phase related to a sym- 
metry breaking of the mean-field Hamiltonian. These 
unphysical second-order phase transitions that appear in 
the slave-boson approach are artifacts of the theory, as 
already observed by Coleman |36| . while the present X- 
boson treatment prevents those spurious phase transi- 
tions. Therefore, we were motivated to study the su- 
perconducting phase diagram systematically by varying 
the hybridization parameter V and our results provide 
three distinct regimes for the superconducting phase di- 
agram as can be seen in Fig. ^ ( a ) f° r small values 
of V the superconductivity shows up in a broad inter- 
val of occupations and the function T c (N t ) have a single 



maximum for our numerical data; (b) as V is increased, 
V ~ 0.2 — 0.4, the range of occupations where there is 
superconductivity is narrowed and shifted to higher val- 
ues. The function T c (N t ) presents two local maxima and 
between them the superconductivity is suppressed, which 
is related to the appearance of a hybridization gap in the 
Pf(oj) and the system becomes an insulator in this re- 
gion; (c) as V increases even more, V ~ 0.4 — 0.5, the 
superconductivity shows up in a smaller range of occupa- 
tion and the superconducting region related to the first 
local maximum is suppressed. As the total occupation 
is raised, the system crosses an insulating region until it 
becomes superconducting. For larger values of the hy- 
bridization parameter (V > 0.5) the superconductivity 
is completely suppressed. 

According to our data for the superconducting phase 
diagram in the first regime, the function T c (N t ) has a 
single maximum, as can be seen in Fig. ^a. For this 
low V regime the /-band density of states is located in a 
narrow region of the bandwidth, since the electrons are 
more localized, the two peaks of the density of states 
are sharp and the distance between them is small, about 
0.02Z3, as can be seen in Figs. 0]a-0]c. Notice that in 
this small interval defined by the distance between peaks, 
the densities of states for both bands are zero and when 
the chemical potential lies in this hybridization gap the 
system presents an insulating behavior unless the con- 
ducting electrons have enough kinetic energy to tunnel 
through the gap. In this low hybridization regime the 
superconducting critical temperature is larger than the 
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FIG. 5: The chemical potential p, as a function of N t . Other 
parameters are the same as in Fig. Q 



distance between the chemical potential p and the bor- 
der band peaks, and since the average kinetic energy of 
the conducting electrons should be proportional to T c for 
a given occupation in the phase diagram, the conduct- 
ing electrons have enough energy to tunnel through the 
hybridization gap and the system does not present an in- 
sulating behavior. Notice that T c is larger for V = 0.1 
because the existence of the superconductivity phase is 
related to a high value for the /-band density of states 
at the chemical potential. 

The phase diagram for the second regime is presented 
in Fig. ^b. The critical temperature T c is non-null only 
when the /-band density of states at the chemical po- 
tential is sufficiently high, what also indicates the range 
of occupations where there is superconductivity. Notice 
the two local maxima for T c and the suppression of su- 
perconductivity between them when p crosses the hy- 
bridization gap. In this region the Matsubara's summa- 
tion defined by Eq. 1(571) is close to zero and, hence, 
there is no non-null solution for T c . Indeed, in Figs. 
0Jd - 01f we present the results for the chemical poten- 
tial localization, the /-band density of states and the 
conducting electrons density of states p c (w) correspond- 
ing to different values of N t . Moreover, Fig. [3] presents 
the behavior of the chemical potential as a function of 
Nt for the three regimes found. Notice that there is an 
abrupt variation of p when it crosses the hybridization 
gap, the slope is larger as the distance between peaks of 
the DOS is larger. Indeed, Fig. 0] shows the appearance 
of a hybridization gap when the chemical potential lies in 
a region between bands where superconductivity is sup- 



pressed because the system becomes an insulator. An 
analogous transition was found by Romano studying 
the same model taking into account s, p, d anisotropic 
pairing and for U — 1, in units of D for the value of 
Nt close to the half filling. Indeed, notice that when the 
values of the hybridization parameter are raised the two 
peaks of the density of states are lowered and the dis- 
tance between peaks is broadened, evidencing two main 
factors responsible for the appearance of an insulator re- 
gion under the constraints of the furnished parameters in 
the X-boson approach: on one hand, when V increases, 
the /-band density takes up smaller values, what causes 
the overall superconducting critical temperature to di- 
minish, since T c is non-null only when the value of Pf{p) 
is sufficiently large, indicating that the average kinetic 
energy of the fermions in the system is decreased; on the 
other hand, since the distance between peaks are broad- 
ened, the demand of kinetic energy required by the con- 
ducting electrons to tunnel through the hybridization gap 
becomes larger. Therefore, as the total occupation is in- 
creased the chemical potential roughly also increases, as 
seen in Fig. 01 and when it crosses the hybridization gap, 
the conducting electrons do not have the required energy 
to tunnel through the gap and the system becomes an in- 
sulator. Again, this result would not be obtained by the 
mean-field slave boson method |15| . since the approach 
breaks down in the Kondo region when Nf — > 1 and the 
upper band cannot be reached. 

For the third regime, superconductivity is constrained 
to a small range of occupation, as can be seen in Fig. 
^c. This regime is characterized by high values of the 
hybridization parameter and the superconducting region, 
related to the first local maximum of T c , is suppressed. 
Therefore, as the total occupation is raised, or charge 
carries are added to the /-band, the system suffers a 
insulator-superconductor transition. Again, our results 
for the superconducting phase diagram can be inferred 
from the behavior of the chemical potential compared to 
several plots of the the /-band density of states. Indeed, 
as can be seen in Figs. Elg-BJi, the /-band DOS is asym- 
metric and the first peak, which is related to the first 
local maximum of T c , is smaller; hence, for this range of 
occupations the values of Pf(p) are not sufficiently large 
and the superconductivity is suppressed. 

It is interesting to compare the results obtained by 
the X-boson approach with those obtained by the slave 
boson approach (see Appendix |A"|) : first note that while 
the effective superconducting interaction W is mapped 
according to W — > D^W for the X-boson method, the 
effective W for the slave boson method is W — > z 2 W, see 
Appendix^ Note that D a is a bounded quantity varying 
from 1 to 0.5 as the /-band occupation raises, while z — > 
as Nf — > 1 and the main effect of the z 2 W interaction 
is to cause a general reduction of T c for the occupations 
where the superconductivity exists. Hence, only for a 
large value of the W parameter we could find a value 
for T c in the slave boson method which is comparable 
to those obtained by the X-boson approach, what means 
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FIG. 6: (a) Slave boson's result for T c as a function of N t . (b) 
2 as a function of Nt. (c) The chemical potential fj, and the 
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W = — 0.5D, and other parameters are the same as in Fig. 



that the slave boson method requires greater values of W 
to produce the same critical temperatures reached by the 
X-boson method, see |3|. Furthermore, Fig. H^a presents 
the superconducting phase diagram for T c as a function 
of N t for V = 0.2 and W = -0.5. Note that as N t 
increases, the chemical potential also increases, but z, 
which is essentially the expectation value of the /-holes, 
tends to zero (Fig. |SJb) causing T c to be null before 
p > E f (Fig. 0c). 

Finally, in Fig. 0we show the parameter g defined by 
the Eq. ^ as a function of N for the three regimes. As can 
be seen in Fig. according to the criteria established in 
0, [24| , as Nt increases we recover the three characteristic 
regimes of the PAM: Kondo, IV and magnetic; the same 
cannot be found for the slave-boson treatment, since it 
breaks down when Nf — > 1. Also notice that the criteria 
is consistent to the picture that as p approaches to Ef 
within an interval of width A a oc V 2 , there exists the 
possibility of charge fluctuations between the bands. For 
V = 0.5, the superconductivity only exists for a range 
of N where the system presents a IV behavior (g > 1). 
As the value of the hybridization is lowered, we prevent 
the possibility of charge fluctuation between bands for 
broader ranges of N t and the superconductivity arises 
even in the regions where the system presents Kondo or 
magnetic behaviors. Since high values of pt(p) usually 
imply in higher values of T C1 the X-boson approach favors 
the superconductivity around the Kondo resonance. Nev- 
ertheless, as fx crosses the hybridization gap, when the 
chemical potential lies in the region between the peaks 
of the density of states, the quantity p c (p) —> 0, and one 
gets that 5 = 0, since it is proportional to the c-band den- 
sity of states at the chemical potential and nothing can be 
inferred about the magnetic order the system. Moreover, 
as Nf — > 1, the parameter g decreases. Indeed, since 
Trkky ~ g 2 and T K oc e~ 1/g [1113 , when g < 1 one 
should expect that the RKKY interaction, which scales 




FIG. 7: Parameter g, as defined by Eq. as a function of 
N t for (a) V = 0.1, (b) V = 0.3, and (c) V = 0.5. Other 
parameters are the same as in Fig. Q 



to ksTfiKKY j to be much larger than IcbTk, suggesting 
that the magnetic ordering at Trkky pre-empts Kondo 
singlet formation, and the system might go from a FL 
paramagnetic state to an RKKY magnetic phase with 
well-localized moments according to the established cri- 
teria. 



VI. CONCLUSIONS 

In conclusion, in this paper we studied within the X- 
boson approach the paramagnetic case of the PAM with a 
phenomenological superconducting term and we obtained 
the phase diagram of the model. Although high values 
of pf{lL) usually imply higher values of T c , what indi- 
cates that the X-boson approach favors superconductiv- 
ity around the Kondo resonance, superconductivity was 
found both for configurations where the system presented 
IV and HF behavior. This behavior is in agreement to 
experimental results, since superconductivity was found 
both in heavy fermion materials as well as in interme- 
diate valence compounds [33, O^]- Moreover, we show 
that higher values of the hybridization parameter implies 
in lower values associated to the maximum of the su- 
perconducting critical temperature, indicating that the 
real charge exchange between bands tends to destabi- 
lizes the Cooper pairing, in agreement to theprevious 
results obtained by the slave-boson approach [j] and by 
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perturbative approach [2(. Furthermore, as the total oc- 
cupation is raised and for larger V, the system presents 
a superconductor-insulator transition, which is related to 
the appearance of a hybridization gap that cannot be ob- 
tained by the slave-boson method, which breaks down in 
this region. These results could be experimentally tested 
since the hybridization coupling V can be increased by 
an applied external pressure. 

Also, the detailed discussion of the magnetic solutions 
of the model within the X-boson approach, the possibility 
of coexistence between superconductivity and magnetic 
order and even how the different symmetries of the super- 
conducting order parameter can alter the results obtained 
in the present paper will be subject of investigation in fu- 
ture works. 
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APPENDIX A: SLAVE BOSON APPROACH 

We apply the Coleman's "slave boson" method 01 
described in Section [Q] to the Hubbard operators in 
the original Hamiltonian given in Eq. [3 plus the term 
Ekk'^k'^k', which is responsible for the Cooper 
pair formation in the singlet state. Here we only as- 
sume a constant hybridization V, as well as site in- 
dependent local energies = Ef i(T . Also, we are 
constrained to the subspace where the identity / = 
b^b + J2a- ftfa i s preserved. Therefore, the operator 



Q = Kyz + J^kcr /k cr/kCT ~ lj is added to the mean-field 
Hamiltonian rewritten in the slave-boson's representation 
and the parameters z and A are determined by minimiza- 
tion of the free energy. In terms of this new representa- 
tion the model Hamiltonian becomes 

k,<7 k,cr 
k.o- 

+J2 z2w w b i b *' 

k,k' 

+A(z-1), (Al) 

where the localized /-energy is renormalized by Ef — 
Ef + A, and the operator b^ is given by b^ = fl ^f_ k j. 
Finally, in a mean-field approximation we find 

(iw n -e k )^(-k,o;„) = V^Q f J(-K^n), (A2) 

(iu! n + £k)^ r c t / , cr (k, LU n ) = -Vy/zFjf^^Uln), (A3) 

(iuj n + E f )F] f>a (k,u> n ) = -Fy^Fj /)CT (k,u/ n ) 

+z 2 A*(k)^(-k,c„), 

(A4) 

(iu n -Ef)gtf(-k,w n ) = VjlQ c J(-\L^ n ) 

+z 2 A k T} fa (k,u n ) 

+ 1 , (A5) 

where the superconducting gap was already defined in 
Eq. |2H1 Also note that the above system of equations 
with the mapping D a = 1, V <-> \fzV and W «-» z 2 W is 
identical to the X-boson's result. 
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